home *** CD-ROM | disk | FTP | other *** search
/ The X-Philes (2nd Revision) / The X-Philes Number 1 (1995).iso / xphiles / hp48hor1 / prob.src < prev    next >
Text File  |  1991-10-19  |  4KB  |  269 lines
  1. %%HP: T(3)A(D)F(.);
  2. @ PROB, by Brian Korver
  3. DIR
  4.   COMBNATION
  5.     \<< \-> N R 'N!/(R!
  6. *(N-R)!)'
  7.     \>>
  8.   BANG
  9.     \<< \-> X 'FACT(X)'
  10.     \>>
  11.   Bin
  12.     \<< 'COMB(n,I)*p^
  13. I*(1-p)^(n-I)'
  14. "Binomial Mass
  15. Function.
  16. Enter n p i"
  17. { ALG V } INPUT
  18. OBJ\-> \-> n p I
  19.       \<< \->STR STR\->
  20. EVAL "P{X=" RCLF
  21. STD SWAP I 4 RND +
  22. SWAP STOF "}=" +
  23. SWAP 4 RND +
  24.       \>>
  25.     \>>
  26.   Bin.D
  27.     \<< '\GS(x=0,I,COMB
  28. (n,x)*p^x*(1-p)^(n-
  29. x))'
  30. "Binomial Distribution
  31. Function.  
  32. Enter n p i"
  33. { ALG V } INPUT
  34. OBJ\-> \-> n p I
  35.       \<< \->STR STR\->
  36. EVAL "P{X\<=" RCLF
  37. STD SWAP I 4 RND +
  38. SWAP STOF "}=" +
  39. SWAP 4 RND +
  40.       \>>
  41.     \>>
  42.   Pois
  43.     \<< 'EXP(-L)*L^I/
  44. I!'
  45. "Poisson Mass
  46. Function.
  47. Enter \Gl (or 'n*p') i"
  48. { ALG V } INPUT
  49. OBJ\-> \-> L I
  50.       \<< \->STR STR\->
  51. \->NUM "P{X=" RCLF
  52. STD SWAP I 4 RND +
  53. SWAP STOF "}=" +
  54. SWAP 4 RND +
  55.       \>>
  56.     \>>
  57.   Pois.D
  58.     \<< '\GS(x=0,I,EXP(
  59. -L)*L^x/x!)'
  60. "Poisson Distribution
  61. Function. (ok \Gl<200?) 
  62. Enter \Gl (or 'n*p') i"
  63. { ALG V } INPUT
  64. OBJ\-> \-> L I
  65.       \<< \->STR STR\->
  66. \->NUM "P{X\<=" RCLF
  67. STD SWAP I 4 RND +
  68. SWAP STOF "}=" +
  69. SWAP 4 RND +
  70.       \>>
  71.     \>>
  72.   Geom
  73.     \<< 'p*(1-p)^(n-1
  74. )'
  75. "Geometric Mass
  76. Function.
  77. Enter n p"
  78. { ALG V } INPUT
  79. OBJ\-> \-> n p
  80.       \<< \->STR STR\->
  81. \->NUM "P{X=" RCLF
  82. STD SWAP n 4 RND +
  83. SWAP STOF "}=" +
  84. SWAP 4 RND +
  85.       \>>
  86.     \>>
  87.   NBin
  88.     \<< 'COMB(n-1,r-1
  89. )*p^r*(1-p)^(n-r)'
  90. "Negative Binomial
  91. Mass Function.
  92. Enter n p r"
  93. { ALG V } INPUT
  94. OBJ\-> \-> n p r
  95.       \<< \->STR STR\->
  96. \->NUM "P{X=" RCLF
  97. STD SWAP n 4 RND +
  98. SWAP STOF "}=" +
  99. SWAP 4 RND +
  100.       \>>
  101.     \>>
  102.   Hyper
  103.     \<< 'COMB(N*p,I)*
  104. COMB(N-N*p,n-I)/
  105. COMB(N,n)'
  106. "Hypergeometric Mass
  107. Function.
  108. Enter N n p i"
  109. { ALG V } INPUT
  110. OBJ\-> \-> N n p I
  111.       \<< \->STR STR\->
  112. \->NUM "P{X=" RCLF
  113. STD SWAP I 4 RND +
  114. SWAP STOF "}=" +
  115. SWAP 4 RND +
  116.       \>>
  117.     \>>
  118.   Normal
  119.     \<<
  120. "Normal Distribution
  121. Unit Transformation.
  122. Enter \Gm \Gs x"
  123. { ALG V } INPUT
  124. OBJ\-> \-> m s x
  125.       \<< "P{X\<=" RCLF
  126. STD SWAP x 4 RND +
  127. "}=\O/(" + '(x-m)/s'
  128. EVAL DUP 3 ROLLD 3
  129. RND + ")=" + SWAP 2
  130. RND 4 FIX 0 1 ROT
  131. UTPN NEG 1 + DUP \->
  132. HLD
  133.         \<< + SWAP
  134. STOF HLD
  135.         \>>
  136.       \>>
  137.     \>>
  138.   UNorm.D
  139.     \<<
  140. "Unit Normal
  141. Distribution.
  142. Enter z"
  143. { ALG V } INPUT
  144. OBJ\-> \-> Z
  145.       \<< "P{Z\<=" RCLF
  146. STD SWAP Z 4 RND +
  147. "}=" + 4 FIX 0 1 '-
  148. RND(Z,2)' \->NUM UTPN
  149. + SWAP STOF
  150.       \>>
  151.     \>>
  152.   DeMoivre.Laplace
  153.     \<<
  154. "DeMoivre-Laplace
  155. Distribution.
  156. Enter n p a b"
  157. { ALG V } INPUT
  158. OBJ\-> \-> n p a b
  159.       \<< 'n*p' EVAL
  160. '\v/(n*p*(1-p))' EVAL
  161. \-> m s
  162.         \<< "'\O/(" '(b
  163. +.5-m)/s' EVAL \->NUM
  164. 3 RND + ")-\O/(" + '(
  165. a-.5-m)/s' EVAL
  166. \->NUM 3 RND + ")=" +
  167. 0 1 '(b+.5-m)/s'
  168. \->NUM 2 RND UTPN NEG
  169. 1 + 0 1 '(a-.5-m)/s
  170. ' \->NUM 2 RND UTPN
  171. NEG 1 + - + "'" +
  172. OBJ\->
  173.         \>>
  174.       \>>
  175.     \>>
  176.   Exp.D
  177.     \<< '1-EXP(-L*I)'
  178. "Exponential Distribution
  179. Function.  (\Gl=1/\Gm)
  180. Enter \Gl i"
  181. { ALG V } INPUT
  182. OBJ\-> \-> L I
  183.       \<< \->STR STR\->
  184. \->NUM "P{X<" RCLF
  185. STD SWAP I 4 RND +
  186. SWAP STOF "}=" +
  187. SWAP 4 RND +
  188.       \>>
  189.     \>>
  190.   Gam.D
  191.     \<< '1-EXP(-L*I)*
  192. \GS(K=0,T-1,INV(K!)*(
  193. L*I)^K)'
  194. "Gamma Distribution
  195. Function.
  196. Enter \Gl t i"
  197. { ALG V } INPUT
  198. OBJ\-> \-> L T I
  199.       \<< \->STR STR\->
  200. \->NUM "P{X<" RCLF
  201. STD SWAP I 4 RND +
  202. SWAP STOF "}=" +
  203. SWAP 4 RND +
  204.       \>>
  205.     \>>
  206.   Weibull
  207.     \<< 'B/A*((X-V)/A
  208. )^(B-1)*EXP(-((X-V)/A)^B)'
  209. "Weibull Density
  210. Function.
  211. Enter \Ga \Gb v i"
  212. { ALG V } INPUT
  213. OBJ\-> \-> A B V I
  214.       \<< \->STR STR\->
  215. \->NUM "P{X=" RCLF
  216. STD SWAP I 4 RND +
  217. SWAP STOF "}=" +
  218. SWAP 4 RND +
  219.       \>>
  220.     \>>
  221.   Weib.D
  222.     \<< '1-EXP(-((I-V
  223. )/A)^B)'
  224. "Weibull Distribution
  225. Function.
  226. Enter \Ga \Gb v i"
  227. { ALG V } INPUT
  228. OBJ\-> \-> A B V I
  229.       \<< \->STR STR\->
  230. \->NUM "P{X<" RCLF
  231. STD SWAP I 4 RND +
  232. SWAP STOF "}=" +
  233. SWAP 4 RND +
  234.       \>>
  235.     \>>
  236.   Cauc.D
  237.     \<< '.5+INV(\pi)*
  238. ATAN(I-T)'
  239. "Cauchy Distribution
  240. Function.
  241. Enter \Gh i"
  242. { ALG V } INPUT
  243. OBJ\-> \-> T I
  244.       \<< \->STR STR\->
  245. \->NUM "P{X<" RCLF
  246. STD SWAP I 4 RND +
  247. SWAP STOF "}=" +
  248. SWAP 4 RND +
  249.       \>>
  250.     \>>
  251.   Beta
  252.     \<< 'FACT(A+B-1)/
  253. (FACT(A-1)*FACT(B-1
  254. ))*I^(A-1)*(1-I)^(B
  255. -1)'
  256. "Beta Density
  257. Function.
  258. Enter a b i"
  259. { ALG V } INPUT
  260. OBJ\-> \-> A B I
  261.       \<< \->STR STR\->
  262. \->NUM "P{X=" RCLF
  263. STD SWAP I 4 RND +
  264. SWAP STOF "}=" +
  265. SWAP 4 RND +
  266.       \>>
  267.     \>>
  268. END
  269.